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We derive the node structure of the radial functions which are solutions of the Dirac equation 
with scalar S and vector V confining central potentials, in the conditions of exact spin or pseudospin 
symmetry, i.e., when one has V = ±S + C, where C is a constant. We show that the node structure 
, for exact spin symmetry is the same as the one for central potentials which go to zero at infinity 

but for exact pseudospin symmetry the structure is reversed. We obtain the important result that 
' it is possible to have positive energy bound solutions in exact pseudospin symmetry conditions for 

^SJ . confining potentials of any shape, including naturally those used in hadron physics, from nuclear to 

_ ■ quark models. Since this does not happen for potentials going to zero at large distances, used in 

' nuclear relativistic mean-field potentials or in the atomic nucleus, this shows the decisive importance 

of the asymptotic behavior of the scalar and vector central potentials on the onset of pseudospin 
symmetry and on the node structure of the radial functions. Finally, we show that these results are 
OO ' still valid for negative energy bound solutions for anti-fermions. 

PACS numbers: 21.10.Hw, 21.60.Cs, 03.65.Pm 



P3 ■ I. INTRODUCTION 

Spin and pseudospin symmetries are SU{2) symmetries of a Dirac Hamiltonian with vector and scalar potentials. 
^ They are realized when the difference, A — V — S, or the sum, T, — V -'r S , are constants. The near realization 
C of these symmetries may explain degeneracies in some heavy meson spectra (spin symmetry) or in single-particle 

energy levels in nuclei (pseudospin symmetry), when these physical systems are described by relativistic mean- field 
■ theories (RMF) with scalar and vector potentials [11, ^ . Recently it was found that nuclear resonant states exhibit 
. I similar features as bound states, namely that in conditions of pseudospin symmetry the same pseudospin quantum 
numbers will be conserved and the pseudospin doublets would have the same energy and width f3|. The spin and 
pseudospin symmetries for a Dirac equation with central Coulomb potentials, together with the node structure of its 
, radial function solutions were also recently discussed in ref. Q. When these symmetries are realized, they decouple 
. . ' the upper and lower components of Dirac equation so its solutions behave, as far as the energy spectrum is concerned, 
J>" I as spin zero solutions of the Klein-Gordon equation with the same vector and scalar potentials Q- For systems whose 
potentials go to zero at infinity, pseudospin symmetry cannot be realized for positive energy solutions but only for 
negative energy solutions ^,3r^\- The reverse is true for spin symmetry. However, for harmonic oscillator potentials, 
. also used as nuclear mean-fields, one is able to find bound solutions when either spin or pseudopin symmetries are 
realized flol - [l2j |. In this paper we will derive the node structure of radial functions for central vector and scalar 
potentials which are finite at the origin and go to infinity as r ^ oo when these potentials satisfy spin and pseudospin 
conditions, independently of their shape. We show in a very general way that for these confining potentials it is possible 
to have positive energy bound solutions for exact pseudospin symmetry, contrary to what happens for potentials going 
to zero at large distances, as is the case of nuclear RMF. This quite general finding means that for these potentials it 
is possible to realize exactly this symmetry in nature, which can be relevant also in particle physics where confining 
potentials like the Cornell potentials are of great interest. Finally, we also show that these results are still valid for 
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bound states of anti-fermions, i.e., that those states can exist in the exact spin symmetry conditions. 



II. NODE STRUCTURE OF THE RADIAL FUNCTIONS FOR CONFINING POTENTIALS 

With no significant loss of generahty, we wih set A = for exact spin symmetry and E = for exact spin symmetry. 
The derivations in this section foUow closely the procedure of Leviatan and Ginocchio [l^ . 
The Dirac Hamiltonian with scalar S and vector V potentials reads 

H = a-pc + /3{mc^ + S) + V , (1) 

where a and /? are the Dirac matrices in the usual representation 

where a are the Pauli matrices and I is the 2x2 unit matrix. The Hamiltonian ([T]) can be written in terms of the 
sum and difference potentials Y, — V + S and A = y — 5 as 

H = a • pc + ;3mc2 + i(/ + /3)S + i(/ - /3)A , (3) 

The general solution of the time-independent Dirac equation Hip = Eip for central potentials is 



4,= 



r 

Mr) 



(4) 



where 



£ is the upper component orbital angular momentum and k — ~k. The angular functions (pKm-i^T^) are the spinor 
spherical harmonics and ^^(r), fkir) are the radial wave functions for the upper and lower components of the Dirac 
spinor respectively. The orbital and total angular momenta can be obtained from k hy £ = \k\ + 1/2(k/|k| — l) and 
j = \k\ — 1/2 . The radial functions ^^(r) and /^(t') satisfy the coupled first-order differential equations 

'^9^ . J^^A{r)h (6a) 
r 



dr 

dr 



df-^ ^fl = ^B{r)g^, (6b) 
r 



where 

A(r) = ^\E + mc^ ^ A(r)] (7a) 
nc 

B(r) ^ ^[E-mc^ - I](r)] . (7b) 
nc 

If the potentials A(r) and S(r) are finite at r = 0, then near the origin one must have g^i^) oc r" and /«('') oc 
where a and (3 are integers greater than or equal to 1. From eqs. ^ one finds, when r — 0, 

f 5„(r) oc r'^+i 

S r , , K > (8a) 

/s(r) cx r'' 



fair) oc r 



-<0. (8b) 
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Setting = r'^g^ and Ta = r'^fa = r '^f^ eqs. © can be written as 



^=r^«A(r)J-^ (9a) 
ar 

'£l^^r-^-B{r)g^ (9b) 
ar 

The functions and J-f^ have the same nodes as g^, and /jj for r > 0. As remarked in ref. [l3| . eqs. ([S]) imply that 
the nodes of and ahernate, i.e., one function radial node is located between two consecutive nodes of the other 
function. One the other hand, if ri and r2 are the the nodes of J-g and C/k respectively, one has 



r—ri 



-A{r^)B{r^)g^[ri) Ff,{r^) ^ Q (10a) 
^A{r2)B{r2)Mr2) ^^(rs) = (10b) 



From and (jlOp and J^ij have an extremum when the other function has a node. Moreover, since we will be 
looking for bound solutions, and must go to zero when r — ^ oo and so these extrema must be a maximum when 
the function is positive and a minimum when the function is negative. Then, from (jlOp . one must have A{r)B{r) > 
at each node. When there is spin symmetry (A — 0), A{r) > 0, so that condition means that one must have 
B{r) = E — mc^ — E > at the nodes of and Ff^ which is to say that the kinetic energy is positive, since S acts 
as a binding potential. This means that in this case all nodes occur within the classically allowed region, a situation 
similar to the case when binding potentials go to zero, as remarked in ref. [13]. In pseudospin symmetry conditions 
(E = 0), and assuming that the confining potentials are positive, one has E > mc^ and therefore B{r) > 0, so that 
A{r) also has to be positive at the nodes of and J-a as a necessary condition for the radial functions to go to zero 
at infinity. Although A cannot be considered strictly a binding potential, A(r) can be related to a position dependent 
effective mass (see, for instance, [31) and so the condition E > —mc^ + A or that A{r) must be positive at the nodes 
is equivalent to the "classical" condition that the effective mass must be positive at those nodes. 

We turn now to the asymptotic behavior of radial functions considering that the potentials A or E are confining 
potentials, i.e., go to infinity as r — >■ oo. To have bound solutions, the radial functions must go to zero at infinity and, 
due to the symmetry of eqs. (O we expect that both ^^(r) and /s(r) will have similar asymptotic behavior so that, 
when r — oo, those equations may be written as 



dg^ 
dr 
dh 



A{r)h (11a) 



-B{r)g, (lib) 

and the corresponding 2nd-order equations as 
d^g^ _ 1 dA{r) dg^ 



dr^ A{r) dr dr 
d^h _ 1 dB{r) df, 



- A(r)B(r) g, ^ -A(r)B(r) g, (12a) 
A{r)B{r) h - ~A{r)B{r) h (12b) 



dr'^ B{r) dr dr 

since the derivative term will go faster to zero than the function term. One sees immediately that the product A{r)B{r) 
must be negative asymptotically in order that the radial functions go to zero at infinity because their second derivative 
must has the same sign as the function. 

One may assume the following asymptotic form for the radial functions (A is a positive constant) 

g.{r) - h{r) ^ e-^f^^^ (13) 

where /(r) is a increasing function of r. Inserting this ansatz into eqs. (1121) we can determine /(r) and the sign of 
the potentials. For instance, if 

A(r) > Cr" (14) 

E(r) > Dr" (15) 
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with a > and C and D are constants, one would have /(r) = r°+^ if both C and D are different from zero or 
/(r) = r°/2+i if either C or D is zero, as would be the case for spin or pseudospin symmetry conditions, respectively. 
For instance, if a = 2 (harmonic oscillator potentials), in spin or pseudospin symmetry conditions, the radial functions 
would behave as Gaussians, as indeed is the case (see [llj). As mentioned before, at the same time we get the sign of 
the potentials from eqs. (jl2p . so that, again in spin or pseudospin symmetry conditions, one has respectively either 
D > (from ([7]) B{r) would be a positive constant, because E > mc^) or C > {A{r) would be a positive constant). 
This is, of course, what one would expect for confining potentials for positive energy solutions, either acting as binding 
potentials like S or as an effective mass like A. 

To obtain a relation between the nodes of g,^ and it is useful to consider the behavior of the product gnfa- From 
eqs. © one obtains 



d{gK.fk) 
dr 



A{r)fl-B{r)gl 



(16) 



Let us consider separately the cases for exact spin and pseudospin symmetries. 



A. Spin symmetry 



At r — > 0, the behavior of the radial functions ([8]) and eqs. ([7]) imply that 
d{gKfk) 



dr 
dr 



-B{Q)gl < K < 
A(0)/| > K>Q. 



On the other hand, when r oo, because Air) is constant for exact spin symmetry, we get 



dr 



'B{r)gi > 



(17a) 
(17b) 

(18) 



Since gnfk is zero at r = 0, eqs. (I17|) give its sign near the origin and the same happens with eq. (jl8p at large r. We 
have 



r ^ gjk<0 {k< 0) 
r^Q gjk > (k > 0) 
r oo < 

Moreover, from ([T6| one has at the nodes ri > and r2 > of and respectively 



dig^fk) 



dr 
d{gK.fk) 



dr 



^B{T^)gl < 



= Air2)fi > 



(19a) 
(19b) 
(19c) 



(20a) 
(20b) 



since, as referred previously, both A{r) and B{r) are positive at the nodes of the radial functions. 

From (jl9p we see that (^k/k changes sign when going from zero to infinity if k > while it is negative near the 
origin as well as at large r if k < 0. This means that g^fk has an odd number of nodes (not counting the origin) when 
K > and an even number of nodes when k < 0. On the other hand, from ([T^ and (^0)) we see that the last node is 
a node of fa because the derivative at the last node must have the opposite sign of the asymptotic derivative. By the 
same argument, from PT|) and d^H]) . when k > the first node (r > 0) must be a fa node, while when k < the first 
node is a g^^ node. Therefore, and since the nodes of fa and g^ alternate, when k > fa must have one more node 
than g^ and when k < the radial functions have the same number of nodes. In summary, denoting by and rig 
the number of nodes of fa and g^ respectively, we have 



K < 
K > . 



(21) 



This relation is known to hold also for central potentials that go to zero at infinity as the nuclear mean-field potentials 
both for fermions and anti-fermions 0, 0, [l3j . 
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B. Pseudospin symmetry 

In pseudospin conditions, the behavior of gnfa at the origin is the same as for spin symmetry, but when r — > oo we 
have 

^^^^ - A(r)/2 < (22) 
dr 

so that now we have 

(23a) 
(23b) 
(23c) 

Equations (j20l) still hold in this case, so using the same reasoning as in the spin symmetry case, we find that now 
the last node must be a node, while the first node will still be a node if k > and a node if k < 0. So in this 
case we will have 

nf^n„ — l K < 

(24) 

Uf — Hg K > U . 

In ref. [ll| was shown that indeed one has such radial node structure for the harmonic oscillator potentials S and 
A in spin and pseudospin symmetry conditions. This structure can be understood by noting that in pseudospin 
symmetry conditions the radial functions interchange roles relative to spin symmetry conditions, so we have (7^ fk 
and K o K = —K. 





9k Ik 


< 


{n 


<0) 






> 


{n 


>0) 


r — > cxD 


.9 k /k 


> 







III. EXISTENCE OF BOUND SOLUTIONS FOR EXACT PSEUDOSPIN SYMMETRY 



A well-known fact about pseudospin symmetry is that it cannot be exact in relativistic mean-field quantum systems 
because a S potential going to zero at infinity cannot not be zero everywhere to provide some binding, so that one 
has E < mc^. As remarked by Leviatan and Ginocchio [l^ the previous statement can be rephrased by saying that 
in order to have bound solutions of the Dirac equation with scalar and vector central potentials of that kind, the 
derivative of g^/s at large r must change sign when r decreases to be able to have a node or at least to go to zero 
at r = 0. For nuclear mean-field potentials this means that, from eq. and since A{r) > everywhere B{r) 
must be able to change its sign from the negative asymptotic value, i.e., one must have B{r) > for some r values, 
or, as noted before, there should be some region in which the kinetic energy is positive, that is, a classically allowed 
region. 

However, if one has confining potentials, one has at hand another mechanism to bind particles, namely by letting 
the effective mass go to infinity as r — >■ 00, which, as referred before, is equivalent to have a confining A potential. In 
this case, one is allowed to let E be zero and still have bound states. This can be seen from eqs. ([T6| and (|22|) since 
a large negative A(r) at large r will become positive when r < rc (where Tc is such that A(rc) = E + rac^) so it is 
clear that there exist values of r < rc such that A{r)f1 > B{r)g-^ — [E ~ mc?)g1^, thus allowing g^fk to change sign 
and fulfilling the necessary conditions to have bound state solutions. 

In refs. [lol [ll| was shown that indeed one can have bound solutions for a harmonic oscillator A potential when 
there is pseudospin symmetry. In this section we have proved quite generally that this is possible for any confining 
scalar and vector central potentials, i.e., that go to +00 when r — > 00. 



IV. BOUND SOLUTIONS FOR ANTI-FERMIONS 



If one considers the spin and pseudospin symmetries for anti-fermions, one has to consider the charge-conjugated 
Dirac equation with central confining scalar and vector potentials. The effect of charge conjugation is basically to 
perform the transformations E ~E, S — s- — A and A — — S, 0|- One gets 

^ + ^h=A{r)h (25a) 
ar r 

^-K^ = -B{r)g, , (25b) 
ar r 
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where 



B{r) 



A{r) 



1 

he 
1 

he 



E + mc^ + T.{r)] 



E - mc^ + A{r)] 



(26b) 



(26a) 



and 5k /k are the radial function for the bound anti-fermions. One can see immediately that the role of the E and 
A potentials is reversed relative to positive energy case, so that now A is the binding potential and E is the "mass" 
potential. In our work concerning the harmonic oscillator potential in the 1 + 1 Dirac equation [l^ . but whose 
conclusions are valid in 3 + 1 dimensions, we have shown that one can have bound solutions in spin and pseudospin 
conditions for anti-fermion for negative oscillator harmonic potentials, i.e., which go to — oo when r ^ oo and that 
the negative energy solutions are such that E < —mc?. Therefore, considering now the case of a general confining 
potential going to — cx), we see that all the reasoning in the previous sections can be repeated with A and E being 
finite potentials at the origin and going to — oo at large distances. One has —E + ttk? > 0, —E — mc^ > 0, and 
thus A{r) and B{r) in (pS)) are positive at the nodes of cjn and when they are constants in pseudospin or spin 
symmetry conditions their asymptotic values are such that their product is negative. Therefore, the node structure 
of 5k fa is the same as ^k fa given by ([2T|) and ([24|. but now in reverse conditions, i.e., in pseudospin and spin 
symmetry conditions respectively. This is consistent with the finding that bound states of anti-nucleons obtained by 
charge conjugation for nuclear mean field potentials have the same radial node structure than their positive energy 
counterparts [6,, . 



In this paper we have derived the node structure of the radial functions of the upper and lower components of a 
spinor which is a bound solution of the Dirac equation with scalar S and vector V potentials such that V = S (spin 
symmetry) and V = —S (pseudospin symmetry), when those potentials are finite at the origin and go to +00 when 
r — >■ 00, independently of their shape. It was shown that the node structure when pseudospin symmetry exists is 
different from the usual node structure in systems of scalar and vector potentials with are finite at the origin and go 
to zero at infinity. However, in spin symmetry conditions, we found that the node structure for confining potentials 
is the same as for potentials which go to zero at infinity. 

We have also shown that, contrary to what happens with potentials which go to zero at large distances, confining 
potentials allow to have bound states in pseudospin symmetry conditions for positive energy states and in spin 
symmetry conditions for negative energy states. 

We believe that the proof presented here of the existence of positive energy bound states for any radial confining 
potential in the limit of pseudospin symmetry is quite relevant, since it uncovers a whole class of potentials which allow 
to realize exactly this symmetry in nature. While it can realized only approximately in physical systems governed by 
interactions going to zero at large distances, as in relativistic mean-field theories for nuclei or in atoms, one can find 
other systems where either harmonic oscillator or linear confining potentials are crucial for explaining their properties, 
as, for example, the Cornell potentials in particle physics. Very recently, solutions of the Dirac equations with scalar, 
vector and tensor Cornell radial potentials in the spin and pseudospin symmetry limit have been obtained, and it was 
shown that the solution in each symmetry limit case can be related by a chiral transformation [Tsj . as was already 
been shown for harmonic oscillator potentials jT^. In that work it is explicitly shown that the radial Dirac equation 
for this problem can be mapped into a Schrodinger-like equation with a harmonic oscillator plus a Cornell potential 
whose solution is presented. This potential is known in particle physics as the Killingbeck potential and was 
also considered recently in the context of the Dirac equation and its solutions in the limit of the spin and pseudospin 
symmetries were obtained [3 ■ Thus we believe that the results presented here can be applied to a wide range of 
physical systems. 



V. 



CONCLUSIONS 
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